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Abstract 

This paper reports the tenth-order QED contribution to the g — 2oi electron and muon from 

^ . two gauge-invariant sets, Set l{g) and Set l{h). In the case of electron g — 2 Set l{g) consists of 

> 

OO ' nine Feynman diagrams which have a fourth-order vacuum-polarization loop containing another 

o ■ 

CN ■ fourth-order vacuum-polarization loop. Set I(/i) consists of 30 Feynman diagrams which have a 

proper sixth-order vacuum-polarization loop containing a second-order vacuum-polarization loop. 
00 ' The results of numerical integration, including mass-dependent terms containing one closed loop 

o 

; of muon, are 0.028 597 (4) (a/vr)^ for Set l[g) and 0.001 685 (13) {a/Txf for Set I(/i), respectively 

^ ' We also report the contributions of Set 1(17) and Set I(/i) to the muon anomaly. Diagrams included 

H ; 

5t , are those containing electron, muon, and tau-lepton loops. Their sums are 2.640 9 (4)(a/7r)^ and 

-0.564 8 (ll)(a/7^)^ respectively The sum of contributions of Sets l{g) and I(/i) containing only 
electron loops are in fair agreement with the recently obtained asymptotic analytic results. 

PACS numbers: 13.40.Em, 14.60.Cd, 12.20.Ds 
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I. INTRODUCTION 



The anomalous magnetic moment g—2 of the electron has played the central role in testing 
;he validity of QED. To match the precision of the recent measurement of the electron g—2 
ij, the theory must include the radiative correction of up to the eighth-order , the 
hadronic contribution 

mm 

, and the electroweak contribution , Mi within the 
context of the Standard Model. As a matter of fact, the largest theoretical uncertainty 
now comes from the not-yet-calculated tenth-order term. Thus, for a more stringent test of 
QED, it is necessary to know, not a crude estimate made in Ref. jll|, but an actual value 
of the tenth-order term. To meet this challenge we launched several years ago a systematic 
program to evaluate the complete tenth-order term [l2, 1^ . 

The tenth-order QED contribution to the anomalous magnetic moment of an electron 
can be written as 

The contribution to the mass-independent term A^^^^ may be classified into six gauge- 
invariant sets, further divided into 32 gauge-invariant subsets depending on the nature of 
closed lepton loop subdiagrams. Thus far, results of numerical evaluation of 18 gauge- 
invariant subsets, which consist of 964 vertex diagrams, have been published jl2,|l^. Some 
of 18 subsets were also analytically calculated 16| and show good agreement with the numer- 
ical results. Several more gauge-invariant sets have been evaluated and are being prepared 
for publication. 

In this paper we report the contribution to A^-^^^ from two gauge-invariant subsets, called 
Set l{g) and Set l{h). Set l{g) consists of diagrams which contain 4th-order vacuum- 
polarization diagrams U^^^ whose internal photon line has an insertion of another 4th-order 
vacuum-polarization loop. This is denoted as 11*^^'^). (See Fig. [H) Set l{h) consists of dia- 
grams which contain proper 6th-order vacuum-polarization diagrams 11*^^) in whose internal 
photon lines a second-order vacuum-polarization loop is inserted. This is denoted as 11*^^'^). 
(See Fig. H) 

The evaluation of the contributions of Set l{g) and Set l{h) to a^e^^ is straightforward if 
the spectral functions of n*^^''*'' and n*-^'^^ are known, since contributions of these vacuum- 
polarization loops can be regarded as superpositions of the second-order g — 2 in which 
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FIG. 1: Typical diagrams of Set l{g) which contain eighth-order vacuum-polarization subdiagrams 
grams belong to this set. 
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FIG. 2: Typical diagrams of Set I(/i) which contain eighth-order vacuum-polarization diagrams 
n(6'2)(g2). Thirty dia grams belong to this set. 

the virtual photon is replaced by massive vector bosons whose distribution is weighted by 



y applicable to our problem 
17 1, and n(4'2) 0. Instead 



the spectral function. Unfortunately this approach is not ful 
since exact spectral functions are known only for n*^^^ n^"^-' 

of pursuing the photon spectral function, we therefore follow an alternative approach [1 
in which eighth-order vacuum-polarization functions H^'*'^^ and n^^'^^ are constructed from 
Feynman-parametric integrals of H*^^^ and Yi^^^ , which are then inserted in the virtual photon 
line of the Feynman-parametric integral of the second-order anomalous magnetic moment 
M^'^\ This procedure leads to the formula 



M2,P = - f dy{l-y)Ii{q^) 
Jo 



y ; 



(2) 



where Il{q'^) represents the vacuum-polarization function inserted into the photon propagator 
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with the momentum q. 

There are two ways to construct 11*^^'^^ One is to treat the entire eighth-order function 
by direct Feynman parametrization. Another is to insert the exact fourth-order spectral 
function of n^^-* 17| in the Feynman-parametric integral of 11'^^^ We choose the second 



approach in this paper. 



For Set l{h) 
tive approach 



an exact spectral function is not yet known. Thus we follow the alterna- 



191 ] which utilizes the parametric representation of the vacuum-polarization 



function 11'^^) itself. Since the parametric integral of 11'^^) is known 19j, it is easy to obtain 
n(6.2) insertion of the second-order vacuum-polarization function or its spectral function 
in n*^^). For simplicity we choose the latter approach. 

Renormalization of the vacuum-polarization function is carried out in two steps. The first 
step is to remove the UV divergences coming from subdiagrams of the vacuum-polarization 
function by the K-operation [19^, supplemented by the R-subtraction fl^ for some dia- 
grams. Then a UV-finite correction, called residual renormalization, is applied to achieve 
the standard on-the-mass-shell renormalization. 

Diagrams discussed in this paper consist of an open lepton line (/i), a closed lepton line 
(Z2), and another closed lepton line (Z3) inserted in an internal photon line of {I2), where (Zj) 
represents electron (e), muon (m), or tau lepton (t). Thus each diagram is characterized by 
a superscript {IMs)- Residual renormalization constants are also denoted by superscripts 
such as (hh), where I2 is an external lepton line and I3 refers to the inner loop. 

The evaluation of Set l{g) is described in Sec. [Tll The evaluation of Set l{h) is described 
in Sec. IIIII Sec. IIV]is devoted to the summary and discussion of our results. For simplicity 
the factor («/vr)^ is omitted in Sees. [ITlandllnl 



II. SET 1(g) 

In this section we consider Set l{g). The renormalized contribution of Set l{g) to the 
magnetic moment of the lepton li is given by the general formula [l^ 

where the first two terms are finite integrals obtained by the K-operation and the last term 
is the residual renormalization term. 
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A. Set 1(5) (^"'^^ 



Let us first consider the case where li,l2,h are all electron lines, namely {hhh) = (eee): 
«(io) [i(^)(eee)] ^ AM^^l^^,) + 2AM(-:i(p,) - 2ABi:$M%- (4) 

This gives a mass-independent contribution to the electron g — 2. The numerical values of 
AM2 p^|j(.p^-i and AM2 p^|j(.p^j obtained by the Monte Carlo integration routine VEGAS [2^ 
are listed in Table [T] and those of AB^p^ and M2|^p2 listed in Table [Tll Inserting these 
values in Eq. (jl]) we obtain 

a(10)[I(^)(eee)] _ q ggg (g)_ (5) 



B. Set l{g)(^^"^) 

The contribution of Set l{g) to Og, in which the inner vacuum-polarization consists of 
muon loop, is given by 

^(10) [i(^)(ee^)] = AM(5I\p,) + 2AM(55(p,) - 2Ai?(7iMi:;L (6) 

where the numerical values of AM2 p^^^p^^ and AM2 p^j^^p^^ are listed in Table [T] and those 
of AB^yl and M2^"p^2 are hsted in Table [TTl m^/rrie = 206.768 282 3 (52) is from Ref. |ll|. 
Inserting these values in Eq. we obtain 

^(10)[j(^)(eem)] ^ Q^^gg 4 (3) x 10"^ (7) 



C. Set l{g)(^'^^) 

The contribution of Set l{g) to a^, in which the outer vacuum-polarization consists of 
muon loop, is given by 

amgY'-^] = ^M^7:l,P., + 2AM(:pt(p,) - 2A</iM(:p1, (8) 

where the numerical values of AM^'^p^^^p^.^ and AM2'^^'j^|p^^ are listed in Table [T] and those 
of AB^p\ and M^'^p^ are listed in Table [TTl Inserting these values in Eq. ([S]) we obtain 

^(10)[j(^)(eme)] ^ q gy^ 5 ^ ^q-4_ (9) 

Terms for (emm), (eet), etc., are even smaller than those for (eem) and {erne). Of course, 
they are easy to evaluate, if needed. However, they are of no physical significance at present. 
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TABLE I: Contributions of diagrams of Set 1(g) of Fig. 1 to the electron g—2 with (eee), (eem), 
and (erne) and muon g—2 with (mee), {mem), and [mme). np is the number of Feynman diagrams 
represented by the integral. All integrals are evaluated by VEGAS 20|] in double precision. 



Integral 


np 


Value (Error) 


Sampling per 


No. of 






including np 


iteration 


iterations 


A 7i#(eee) 
^^'^2,P4a(P4) 

^^'^2,P4f)(P4) 


3 
6 


0.028 534 (5) 
0.005 798 (2) 


1 X 
1 X 


10^ 
10^ 


50 
50 


^^'^2,P4a(P4) 
^^'^2,P4f)(P4) 


3 
6 


0.171 9 (2) X 10-4 
0.002 1 (0) X 10-4 


1 X 
1 X 


10*5 
10*5 


20 
20 


. (eme) 
^^'^2,P4a(P4) 

^^*^2,P4b(P4) 


3 
6 


0.052 2 (1) X 10-4 
0.044 6 (1) X 10-4 


1 X 
1 X 


10^ 


20 
20 


^^'^2,P4a(P4) 
^^'^2,P4b(P4) 


3 
6 


0.533 54 (16) 
2.219 96 (17) 


1 X 10^ 
1 X 10^ 


1 X 10* 
1 X 10* 


150, 50 
150, 50 


A nA'mem) 
^^*^2,P4a(P4) 

A A/fCnief") 
^^'^2,P4f)(P4) 


3 
6 


0.059 619 (60) 
0.002 157 (12) 


1 X 
1 X 


10^ 
10^ 


20 
20 


^^'^2,P4a(P4) 
A A/f 

^^'^2,P4f)(P4) 


3 
6 


0.139 86 (19) 
0.131 30 (23) 


1 X 
1 X 


10^ 
10^ 


20 
20 



D. Contribution to muon g — 2 from Set l{g) 

The leading contribution of Set l{g) to the muon anomaly comes from the diagrams 
in which both vacuum-polarization loops consist of electrons, namely, 

where the numerical values of AM^^^^^^^p^,^ and AM^^^^^^^p^.^ are listed in Table [T] and those 
of AB^^p^ and 2* are listed in Table [Tll Inserting these values in Eq. (11 01) we obtain 

a(io)[i(^)(™ee)] ^ 2.351 53 (23). (11) 
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TABLE II: Auxiliary integrals for Sets l{g) and I(/i). Some integrals are known exactly. Other 
integrals are obtained by the integration routine VEGAS. The superscript i-*4 stands for the sum 
of the parametric constructions of the fourth-order vacuum-polarization functions P4a and P46. 
For the case of Ai?2,P4, the fourth-order vacuum polarization is inserted into AB2 by using the 
Kallen-Sabri spectral function. The relation between P4 here and the sum of P4a and P46 is given 

by Ai?2,P4 = Ai?2,P4a + Ai?2,P4fe — AB2AB2,P2- 

Integral Value(Error) Integral Value(Error) 



IVI2P2 


0.015 687 421... 
0.005 197 (1) X 10"^ 
1.094 259 6 (0) 




0.076 401 785... 
0.027 526 (1) X 10" 
3.134 97 (14) 


AB2 

ABi% 

^B^2pI 
ABi^i 


0.75 

0.063 399 266... 
0.094 054 (1) X 10"^ 
1.885 766 (77) 


A o(ee) 
^-°2,P4 

. „(em) 
^-°2,P4 

. „(me) 
^-°2,P4 


0.183 666 5 (18) 
0.338 738 (12) x 10 
2.438 91 (23) 


^^'^2,P4{P2) 
^^'^2,P4(P2) 


0.013 120 (1) 

0.078 454 (8) x 10"^ 


^^'^2,P4(P2) 
^^'^2,P4(P2) 


0.047 678 (13) x 10 
1.579 51 (4) 


A o(ee) 
^-°4,P2 

. „(em) 
^-°4,P2 

. „(me) 
^-°4,P2 


-0.314 320 (10) 
-0.915 7 (47) x 10-"^ 
-3.420 4 (72) 


^-^4,P2 

^-^4,P2 

A r ("le) 
^-^4,P2 


0.200 092 (14) 
0.116 1 (64) X 10-"^ 
3.121 3 (63) 



Contributions of terms of type [mme) and {mem) are 

W^) ^"""'^] = 0. 194 64 (29) , (12) 

and 

aaO)^I(^^(mem)j _ q qq^ (61), (13) 

respectively. 

Terms involving tau-Iepton loop are three orders of magnitude smaller than f|TT]) . but 
greater than the uncertainty of (ITT!) . Contributions of terms of type (met), (mte), (mmt), 
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{mtm), and {mtt) are 



respectively. 





\H / J 


= 0.001 236 (2), 






= 0.001 485 (4), 






= 0.000 893 (2), 






= 0.000 663 (1), 






= 0.000 141 (1), 



(14) 
(15) 
(16) 
(17) 
(18) 



III. SET l{h) 

For Set \{h) the formula for the renormalized quantity, including residual renormalization, 
takes different forms depending on whether one follows the original K-operation [19], which 
subtracts only the UV-divergent part of mass- renormalization constant, or the recently de- 
veloped R-subtraction method [14], which subtracts the entire mass-renormalization term. 
We follow here the latter approach which affects the definition of pg'^p2) i = c,d. 
Of course, these two methods are equivalent and give the same results. 



A. Set 

The renormalized contribution of Set l{h) with (eee) to is given by 

ai^°)[i(/.)(-^)] = x:am(;J(P2) 

i=a 

+ 10AB,ABi%M^% - 2iAL[% + ABi%)M^%, (19) 



where AM^'^p^f^p^^ AL^^^p2, and AB^p2 are obtained from the lower-order relations 



AM2,P4 = AM2,P4a + 2AM2,P46, 

AL4 = AL4^ + 2AL4C + AL4i + 2AL4^. 

A54 = AB^a + AB^b, (20) 
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by insertion of a vacuum-polarization loop P2 in each of two photon lines of P^a, P^b, 
and AB4. Substituting numerical values listed in Tables HTl and IIIII in Eq. f[T9|) we obtain 



B. Set 

The contribution of Set l{h) to ag, in which the inner vacuum-polarization loop consists 
of muon, is given by 

a(^°)[i(/.)(-)] = x:^<mP2) 

i=a 

+ 10AB2AB^^;;Im^% - 2{AL[yl + ABtp^)Mi%. (22) 

Numerical values of AM^^^'^p^jpg) listed in Table III, and auxiliary quantities are listed in 
Tables II. Inserting these values in Eq. ( |22|) we obtain 

a{10)[I(/^){eem)| _ _q_233 5 (I3) x 10"^ (23) 

C. Set I(/i)("™") 

The contribution of Set \{h) to Og, in which the outer vacuum-polarization consists of 
muon loop, is given by 

h 

' A A/, 

'2,P6j(P2) 



- AAB^AM^^^^ll^^- AAB'ip.AM^^p, 

+ lOAB^ABiy^M^J^ - 2(ALi7) + AB^^;l)M!,yl (24) 

Numerical values of AM2^^'^|p2) are listed in Table III, and auxiliary quantities are listed in 
Tables II. Inserting these values in Eq. ( I24|) we obtain 



a. 



(10) 



e 



= 0.127 9 (2) X 10-^ (25) 



Terms for {emm), (eet), etc. will be even smaller. Thus, they are of no physical signifi- 
cance at present. Of course it is easy to evaluate them, if needed. 
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TABLE III: Contributions of diagrams of Set 1(h) of Fig. 2 with (eee), (eem), and (eme). nj? is 
the number of Feynman diagrams represented by the integral. All integrals are evaluated in double 
precision. 



Integral 




Value (Error) 


Sampling per 


No. of 






including np 


iteration 


iterations 



A 7M-{eee) 
^'"^2,P6a(P2) 

A7M-(eee) 
^-"^2,P6c(P2) 

^-'^^2,P6d(P2) 

^-'^^2,P6e(P2) 
A T r(eee) 



A/M-(eee) 
^-'^-'2,P6/(P2) 

^-'^-'2,P6fl(P2) 

^-'^^2,P6fe(P2) 

. ^(eem) 
^-"^2,P6a(P2) 

. ^(eem) 
^-"^2,P6fe(P2) 

^-'^^2,P6c(P2) 

^-'^^2,P6rf(P2) 

^-'^-'2,P6e(P2) 

^-'^-'2,P6/(P2) 

^-'^^2,P63(P2) 
A 7K7-(eem) 
^-'^-'2,P6fe(P2) 

^-'^^2,P6a(P2) 

^-'^-'2,P66(P2) 

^-'^^2,P6c(P2) 
A /^^(eme) 
^-'^-'2,P6d(P2) 

^-'^-'2,P6e(P2) 

^-'^^2,P6/(P2) 

^-'^-'2,P6g(P2) 



4 
2 
4 
4 
8 
4 
2 
2 



0.003 269 (2) 
0.002 367 (2) 
-0.003 602 (2) 
-0.006 761 (3) 
0.072 639 (9) 
-0.051 169 (4) 
0.028 956 (4) 
0.003 673 (3) 



4 


0.000 99 (7) X 10-^ 


1 X 10^ 


50 


2 


0.001 08 (5) X lO-'' 


1 X 10^ 


50 


4 


-0.006 60 (5) X 10-"^ 


1 X 10^ 


50 


4 


-0.014 12 (14) X 10-^ 


1 X 10^ 


50 


8 


0.463 55 (84) x 10"^ 


1 X 10^ 


50 


4 


-0.476 36 (57) x 10"^ 


1 X 10^ 


50 


2 


0.140 65 (50) X 10-^ 


1 X 10^ 


50 


2 


-0.198 15 (63) X 10-^ 


1 X 10^ 


50 


4 


0.073 23 (7) X 10-^ 


1 X 10^ 


50 


2 


0.047 15 (5) X lO-'' 


1 X 10^ 


50 


4 


0.015 22 (6) X 10-^ 


1 X 10^ 


50 


4 


-0.061 12 (6) X 10-^ 


1 X 10^ 


50 


8 


0.398 96 (16) X 10"^ 


1 X 10^ 


50 


4 


-0.150 36 (5) X lO-'' 


1 X 10^ 


50 


2 


0.122 27 (5) X 10-^ 


1 X 10^ 


50 


2 


0.018 79 (-1) X 10^1 


1 X lO'^ 


50 



X 10^ 
X 10^ 
X 10^ 
X 10^ 
X 10^ 
X 10^ 
X 10^ 
X 10^ 



150 
150 
150 
150 
150 
150 
150 
150 
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D. Contribution to muon g — 2 from Set I(/i) 



FORTRAN programs for the electron g — 2 can be readily converted to the muon case 
by changing one or two parameters. In the case (mee) the renormalized contribution of Set 
\{h) to is given by 

a?°) [!(/.)(-)] = X:AM(--}^,) 

i=a 

+ lOAEaAE^^;^^^^^^) - 2(ALi52 + ^Bi%)M^jl (26) 
Substituting numerical values listed in Tables [TTl and [IV] in Eq. (126!) we obtain 

a(10)[I(;^)(mee)j _ _q ygg gg ^5g^_ (27) 

Contributions of terms of type {mme) and {mem) are 



and 



a^-'[l{hy"'"''>] = 0.253 70 (57), (28) 



q(10) (mem)j _ _q gg^ (g2), (29) 



respectively. Contributes involving tau lepton is small but not negligible. Contributions of 
terms of type (met), (mte), (mmt), (mtm), (mtt) are 



respectively. 



'"^i)] = -0.001 303 (17), (30) 

™*^)] = 0.003 281 (5), (31) 

mmt)^ = -0.000 612 (4), (32) 

™*™)] = 0.000 746 (2), (33) 

= 0.000 026 (1), (34) 



IV. SUMMARY AND DISCUSSION 

The contribution of Set l{g) to the electron g—2, including mass- depending terms where 
one of lepton loops is a muon loop, is the sum of ([5]), ([7]), and ([9]): 

5 



a(io)[I(^)] = 0.028 592 (6) . (35) 
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TABLE IV: Contributions of diagrams of Set 1(h) of Fig. 2 to with (mee), (mem), and [mme). 
np is the number of Feynman diagrams represented by the integral. All integrals are evaluated in 
double precision. 

Integral np Value (Error) Sampling per No. of 

including np iteration iterations 



. ^(mee) 

^-"^2,P6c(P2) 

. ^(mee) 
^-'^^2,P6d(P2) 

. ^(mee) 
^-'^^2,P6e(P2) 

A 7^f(mee) 
"-'"2,P6/(P2) 

^-'^-'2,P63(P2) 
* ^(mee) 
^-'^^2,P6fe(P2) 



^-"^2,P6a(P2) 

^-'^^2,P6fe(P2) 

^-'^^2,P6c(P2) 
. ^(mem) 
^-'^^2,P6rf(P2) 



^-'^-'2,P6e(P2) 
^-'^-'2,P6/(P2) 
^-'^^2,P63(P2) 
^-'^-'2,P6fe(P2) 



^-'^^2,P6a(P2) 
^-'^-'2,P66(P2) 
^-'^^2,P6c(P2) 
^-'^-'2,P6d(P2) 
^-'^-'2,P6e(P2) 



2,P6/(P2) 

2,P6g(P2) 
^(mme) 
2,P6^(P2) 



4 


4.039 70 


(23) 


1 


X 


109 


200 


2 


2.198 75 


(8) 


1 


X 


109 


200 


4 


1.738 74 


(9) 


1 


X 


109 


200 


4 


-3.303 39 


(8) 


1 


X 


109 


200 


8 


7.478 07 


(29) 


1 


X 


109 


200 


4 


-12.321 11 


(19) 


1 


X 


109 


200 


2 


-1.867 04 


(23) 


1 


X 


109 


200 


2 


6.007 99 


(9) 


1 


X 


109 


200 


4 


0.003 32 


(13) 


1 


X 


10^ 


50 


2 


0.002 236 (77) 


1 


X 


10^ 


50 


4 


0.002 148 (25) 


1 


X 


10^ 


50 


4 


-0.000 639 (80) 


1 


X 


10^ 


50 


8 


0.289 64 


(50) 


1 


X 


10^ 


50 


4 


-0.231 04 


(18) 


1 


X 


10^ 


50 


2 


-0.047 18 


(51) 


1 


X 


10^ 


50 


2 


0.014 71 


(7) 


1 


X 


10^ 


50 


4 


0.206 70 


(19) 


1 


X 


10^ 


50 


2 


0.134 50 


(15) 


1 


X 


10^ 


50 


4 


0.036 45 


(14) 


1 


X 


10^ 


50 


4 


-0.169 33 


(13) 


1 


X 


10^ 


50 


8 


1.077 82 


(44) 


1 


X 


10^ 


50 


4 


-0.451 81 


(13) 


1 


X 


10^ 


50 


2 


0.281 06 


(11) 


1 


X 


10^ 


50 


2 


0.146 48 


(9) 


1 


X 


10^ 


50 
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The contribution of Set l{h) to the electron g—2, including mass-depending terms where 
one of lepton loops is a muon loop, is the sum of fl2T]) . fl25l) . and fl23l) : 

5 



[l{h)] = 0.001 685 (13) (^-j . (36) 



Thus the contributions of Set l{g) and Set l{h) to the electron q—2 are small and within 
computational uncertainties of several other gauge- invariant sets [j]. 

The contribution of Set l{g) to the muon g — 2, including terms where lepton loops are 
electrons, an electron and a muon, an electron and a tau lepton, or muon loops, is the sum 

of (nn), (HI, (m, (iisi), (Hzd, and ©: 

a(^°)[%)] = 2.630 86 (37)(^)', (37) 

noting that Eq. holds for the case {mmm), too. 

The contribution of Set l{h) to the muon g — 2, including terms where lepton loops are 
electrons, an electron and a muon, an electron and a tau lepton, or muon loops, is the sum 
of dlZD, dH, ([29D, dnD, il, dMD, and din): 
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[1(h)] = -0.564 8 (11) , (38) 



noting that Eq. fl2T|) holds for the case [mmm), too. 

Thus the contributions of Sets 1((7) and 1(h) to the muon (7—2 are also very small compared 
with the total contribution 663 (20) (a/vr)^ of 18 gauge-invariant sets, which include all dom- 
inant terms containing light-by-light-scattering subdiagrams and/or vacuum-polarization 



subdiagrams 12(]. Concerning the total contribution to the muon g — 2, see also 211 ]. 



Recently an analytic value of the asymptotic contribution of the sum of Set l{g) and Set 



l{h) to has been obtained: 22| 



a^^^lKg + h)^""^^^ : analytic-asymptotic] = 1.501 12 (^^) . (39) 

This is in fair agreement with our result 

4^°)[I((7 + h)^""''^] = 1.560 70 (64) (^^)' . (40) 

Note, however, that the difference between fl39l) and fl40p is 93 times larger than the estimated 
uncertainty of numerical integration. This may be attributed to the 0{me/m^) term not 
included in the analytic result. 
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